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and to the matter Lagrangian. The model is compared to the Lambda Cold Dark 
Matter concordance model and to the observational data using the "gold" SNela 
sample of Riess et. al. (2004). An excellent fit to the data is achieved. We present a 
toy model designed to demonstrate that such a new, possibly fundamental, principle 
may be responsible for the recent period of cosmological acceleration. Observational 
constraints remain to be imposed on these models. 
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1. Introduction 

One of the most profound discoveries of observational physics is that the universe is 
accelerating in its expansion [H [2j [31 [U, [5l [6] . There have been many attempts to explain 
this late-time acceleration, for example, a pure cosmological constant, dark energy 
associated with some new scalar field and modified gravitational theories, although 
all current models require some level of fine-tuning and none are considered to be a 
complete explanation Lp. Whatever is responsible for the current acceleration may arise 
from some completely new physical principle. This is the possibility we consider in this 
paper. Our goal is to construct a toy model that represents a late-time accelerating 
Universe using a new, possibly fundamental, principle. As our guiding principle, we 
hypothesize the existence of a minimal curvature scale in gravity. 

In a Friedmann, Robert son- Walker (FRW) space-time, without cosmological 
constant A and with only standard matter sources such as dust and radiation, the 
universe will always decelerate as it expands. One way to avoid this is to add matter to 
the system that violates the strong energy condition (SEC). In a cosmological context 
this violation constitutes the addition of matter sources satisfying the equation of state 
p < — l/3p. A second possibility is to explicitly remove flat space-time as a solution to 
the theory. In this case the vacuum of the theory, which is approached at late times 
as the energy density in matter fields becomes more and more dilute, is not Minkowski 
space-time, but instead an accelerating Universe [91 [TOj [11] . To remove flat spacetime 
as a solution we hypothesize the existence of a minimal curvature in our underlying 
fundamental theory. The simplest example of this is, of course, to introduce a bare 
cosmological constant into General Relativity. However, in principle there may exist 
many other ways to achieve this result. Indeed, it appears that many accelerating 
cosmological models derived from modified gravity theories contain such a minimal 
curvature [321 M, E] . 

The idea of a minimal curvature scale in gravity mirrors that of a maximal curvature 
scale. In the literature many authors have considered this possibility and used it 
to remove the curvature singularities of General Relativity by bounding curvature 
invariants from above at the level of the classical action [15]- [26]. In the case of 
singularity removal, it is necessary to bound all curvature invariants in order to cover 
all possible physical situations in which such a singularity may occur. 

By contrast, in the case of a minimal curvature approached at late times in a 
homogeneous, isotropic universe, symmetry implies that it is only necessary to bound 
the Ricci scalar R from below. Hence, unlike in the case of a maximal curvature 
hypothesis, we shall see that one may implement a minimal curvature by using a modified 
Brans-Dicke theory where the Brans-Dicke field couples non-minimally to the matter 
Lagrangian. 

Within this context we demonstrate that the existence of the minimal curvature 
(MC) produces a Universe that evolves from a matter dominated period to an 

| For a recent reviews on the subject see [8] . 
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accelerating phase mimicking the A-Cold-Dark-Matter (ACDM) model. We emphasize 
that the model presented here is only a toy construction of the late Universe. The 
model is not intended to provide a consistent cosmology from the time of Big-Bang 
Nucleosynthesis (BBN) until today. It is unlikely that the precise model presented 
here is compatible with solar system experiments and the tight constraints on the time 
variation of Newton's constant. However, the model does provide an example of how 
the postulated existence of a minimal curvature scale in gravity can provide a new 
mechanism to generate cosmological acceleration of the late Universe. Furthermore, 
the model may capture features of a possibly more fundamental theory that admits a 
minimal curvature scale 0- 

In Section [2j we describe the minimal curvature construction, first by using a toy 
example and then by using a class of modified Brans-Dicke theories. We solve the 
equations of motion for this example and demonstrate how the Universe evolves from 
a matter dominated phase to an accelerating period as the curvature approaches its 
minimal value. In Section [3j we compare the MC model with ACDM and to the 
supernovae (SNela) "gold" sample of [27]. Finally, we comment on the possibility 
of constructing more realistic models that satisfy the limiting curvature hypothesis 
and offer our conclusions and speculations in Section |H In Appendix A, we provide 
a detailed analysis of the vacuum MC theory. In Appendix B, we construct an Einstein 
frame description of the vacuum theory and compare it to the MC vacuum. 

2. The Minimal Curvature Construction 

Our goal is to construct theories in which a certain physical quantity is bounded from 
below. Before leaping directly into our model, it is instructive to consider an example 
of how a similar effect may be achieved in a simpler theory - the bounding of velocities 
from above in Special Relativity by the speed of light [20] . 

The Newtonian action for a free particle of mass m in motion is 



In this classical theory the velocity of the particle is without bound. Now let us implement 
one of the fundamental consequences of Special Relativity: To ensure that the speed of 
this particle is limited by the speed of light we introduce a field cj){t) which couples to 
the quantity in the action that we want to bound (x 2 ) and has a "potential" U{<p). The 
resulting action is 




(1) 




(2) 



The variational equation with respect to <p 
., dU 




(3) 



S An interesting example of a minimal curvature scale occurs in a certain classical limit of quantum 
gravity [131 QJ]. 
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ensures that x is bounded, provided dU /dtp is bounded. Note the absence of a kinetic 
term for tp in the action, and hence, the reason the word potential appears in quotes 
above. In order to obtain the correct Newtonian limit for small x and small we take 
U(4>) proportional to 2 . In the Newtonian limit the action (j2j) reduces to ([T]). A simple 
potential satisfying the above asymptotics is 

™ = T+24, ' « 

Integrating out (f) yields (up to an irrelevant constant) the action for relativistic particle 
motion: 



SsR = m I dt vi — x 2 . (5) 



The above model provides a powerful example of how a toy construction based on a 
fundamental principle - the existence of a universal "speed limit" - can capture features 
of a more fundamental theory. 

We now use a similar construction to model the existence of a minimal curvature 
(MC) scale in gravity. Because we are interested in late time cosmology, we need only be 
concerned with bounding one curvature invariant, the Ricci scalar R. In direct analogy 
with our example from Special Relativity, we introduce a scalar field <p that couples to 
the quantity we wish to bound, R, and a "potential" function V(ip) 

r I ' M 2 \ 
Smc = Afj - 7<pR + V{ip) , (6) 



where M p \ = (87rG) -1 / 2 is the reduced Planck mass and both <p and 7 posses dimensions 
of mass, [y?] = [7] = M. The vacuum theory is equivalent to a Brans-Dicke theory with 
Brans-Dicke parameter ujbd = 0. This is seen by re- writing the action ([6]) in terms of a 
new field 

*=(^-w), (7) 

so that the action becomes 

S BD =jd i x^{®R + V($)) . (8) 

Ordinarily, such a theory can be re-cast as a purely gravitational theory with 
Lagrangian C = y/—g~f(R) (see, e.g. [29]); however, this is not possible for all forms 
of the potential V. 

For reasons that will become clear shortly, we allow the field (p to couple non- 
minimally to matter. The non-minimal coupling yields the matter action, 

S M = J d 4 x y^g C M (ipu V, fiv/Mpi)^) , (9) 

where C m is the Lagrangian made up of whatever matter fields ipi are in the theory. In 
the case were if) represents a dark matter Dirac spinor, the field p couples non-minimally 
only to the dark matter sector and need not couple to baryons. In this case it is possible 
to avoid constraints on such a coupling from solar-system and table-top tests of gravity. 
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In string theory, non-perturbative string loop effects do not generically lead to universal 
couplings, allowing the possibility that the dilaton decouples more slowly from dark 
matter than ordinary matter (see, e.g. [3QI ED E21 E31 EH EHJ 136]). 

This coupling can be used to address the coincidence problem, since the acceleration 
is triggered by the coupling to matter. For the purposes of our toy construction, we do 
not distinguish between baryonic and dark matter in the remainder of our discussion. 

Note that in the case of f(R) theories which are conformally identical to models 
of quintessence in which matter is coupled to dark energy with a large coupling, 
this strong coupling induces a cosmological evolution radically different from standard 
cosmology [37]. Similar difficulties may arise in the model presented here, however, we 
have yet to investigate this issue. 

The matter stress-energy tensor is given by 

T,^--L^#^. (10) 

We assume a perfect fluid 

T^ u = (p M + Pm)u^u u + P M g^u , (11) 

where u a is the fluid rest-frame four- velocity, and the energy density pu and pressure 
Pm are related by the equation of state Pm = wpM- Because we are focusing on the late 
Universe we shall ignore the presence of radiation and consider only a matter density 
Pm, which redshifts with expansion in the usual manner (with the exception of the 
non-trivial <p dependence) 

pu (12) 

In the above, / is a function describing the non-minimal coupling of the field <p to the 
matter Lagrangian. Such couplings in the context of ordinary scalar-Einstein gravity 
were studied in [31] where it was found that this coupling can be made consistent 
with all known current observations, the tightest constraint coming from estimates of 
the matter density at various redshifts. This coupling plays a critical role in Modified 
Source Gravity, introduced in [38J. 

Variation of the total action S to t = Sue + Sm with respect to the metric tensor, 
bStot/bg^ = yields the modified Einstein equations 

1 

r 

+ aV^V u ip - ag^Uip -V g fll/ = 8nG (13) 

where □ = g^V^V,, and we have introduced a = 2j/M? l . Variation with respect to 
the field <p gives 

dV 1 0/te/Mpi) 
1 R = 7T + 7 a Pm ■ ( 14 ) 

Equation Q14p is the key to imposing the limiting curvature construction. It is clear that 
the curvature R will remain bounded and approach a constant curvature at late times 



[l - a(p)R^ - -(1 - a<p) Rg^v 
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(t — > oo) as long as V + f'pu —* constant; this is the essence of the construction. We 
assume a flat (k = 0) Friedmann-Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t){dr 2 + r 2 dVl 2 } , (15) 

where a(t) is the scale factor of the Universe and dQ is the line-element of the unit 
2— sphere. Defining the Hubble parameter by H = a/ a and substituting the metric 
ansatz ( fl5l) into (fl3l) and (fl4l) gives the generalized Friedmann equation (the 00- 
component of ffl3|) ) and the equation of motion for tp: 

/fa/M pl )p£ } 



ZM^H 2 - Q 1V H 2 - Q-fHip + V(<p) 



foa 3 



6 7 (2tf 2 + H)- V'(cp) = f[J£-\M- l (16) 



where a prime denotes differentiation with respect to (p, /o = f((po/M p \), p and 
are the values of p and pu today. 

By considering the asymptotics of our cosmology at early and late times, we can 
constrain the forms of the functions V(p) and f(<p). We require that the effective 
Newton's constant for our theory remain positive definite so that gravity is always 
attractive. This imposes a constraint on $ > 0, V t in equation (jSJ). There are 
rather strong constrains on the time variation of Newton's constant from the period of 
nucleosynthesis until today (roughly, G/G < 0(1O~ 10 — 10~ 13 yr^ 1 ) [39]). For the time 
being, we will allow ourselves to ignore this constraint in order to produce a toy model 
capable of realizing the MC conjecture. Furthermore, because we are only interested 
in the behavior of the Universe from the matter dominated epoch until today, we have 
ignored the presence of radiation. The absolute earliest our theory is valid is up to the 
period of equal matter and radiation domination t eq . For specificity, by early times we 
refer to times near the time of photon decoupling at a redshift of z ~ 1100, during which 
the Universe is typically already well into the matter dominated regime. 

At late times (and low curvatures) we want to bound the Ricci scalar R from below. 
This will constitute a successful example of a model obeying the minimal curvature 
hypothesis. To bound the curvature we use equation (JHJ). It is clear that the curvature R 
will remain bounded if we bound V' + f'pM, where the prime denotes differentiation with 
respect to (p. Hence, we require r y~ 1 {V' + f Pm) — > 7?. at late times, where we denote the 
hypothesized minimal curvature scale by TZ. We anticipate that, by construction, there 
will be a late-time attractor that is a constant curvature space-time with R = TZ. This 
attractor is not an actual solution to the equations of motion. The above considerations 
restrict the functional forms of potential V and the non-minimal coupling function /. 
Integration singles out a class of theories that must obey V((p) + f{<p)pM ~^ jT^p as 
t —>■ oo. The simplest forms for V and / obeying the above constraint are the linear 
functions 

V%0=MV, /fa) = 17", ( 17 ) 



M 



pi 



The Accelerating Universe and a Limiting Curvature Proposal 



7 



where /i is, in principle, another free parameter with dimensions of mass. However, we 
will take /i — 7, so that 1Z = 7 2 , and 7 is the only free parameter in the theory [jj. We 
now rescale time t = tHo so that today to = 1 and introduce the following dimensionless 
quantities 

H S l VS *i S -l- (18) 
a H M pl M p i 

Q _ Pm _ Pm r _ R hQ x 

^ = ^-3^5'^-^' (19) 

and take ao = 1, where ao is the value of the scale factor today. 
In terms of the dimensionless quantities the EOM become 

^-2^-2^0+^ = 1^, (20) 

3 f a 3 

7 (2tf 2 + H)-YM = f^ ) ^L j (21) 



The successful implementation of the minimal curvature hypothesis is now apparent. 
Recasting equation (1211) in terms of the curvature scalar R = 6(2H 2 + H), and 
substituting in our choices for V and / ([1] 



i? = 7 2 + 3 ^. (22) 

We see that, as the Universe expands and the matter term dilutes, we asymptotically 
approach the minimal value of the curvature R = 1Z = 7 2 . 

It is both interesting and surprising that the solution to Eq. (1221) reduces to the 
simple case of ACDM plus an arbitrary amount of dark radiation which may have either 
positive or negative effective energy density. Most notably, this model arises in Randall- 
Sundrum brane cosmology which has been extensively studied in the literature 0. The 
Friedmann equation derived from the Randall- Sundrum model for a flat Universe is 

R2= 3i% + 36?W| + ^' (23) 

where M 5 is the five dimensional Planck mass and Ca~ A is the so called dark radiation 
term, since it scales like radiation, but it's origins are purely gravitational and it does 
not interact with standard matter [15]. At low energies (when the energy density is 
much less than the critical brane tension), the p 2 term can be safely neglected. The 
main observational restrictions on the dark radiation term come from the acoustic 
scale at recombination (see, e.g. [46]), and from the amount of total growth of density 
perturbations in the non-relativistic matter component from the time of equal matter 
and radiation until the present day [13 HH]. As a result, the density of dark radiation 
cannot be significantly larger than the present CMB energy density □. 

|| It is interesting to note that if we do not take /1 = 7 and then take the limit that 7 — > the theory 
resembles the action of the modified source gravity models studied in [38j . 

For some reviews see, e.g. [43] [44] . 
+ Note, the constraints from BBN are even stronger [49] . 
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Making use of (fTSD and tHM. Eq. 



2 Vl M Vt R 



may be recast (neglecting the p 2 term) as 

(24) 



where we have included a cosmological constant term Q\ = A/SM^H 2 , 
Qc + includes contributions from both the dark and ordinary radiation, 
we find 



H 



-3Q 



M 



2ft, 



2a 3 a 4 
Constructing the Ricci Scalar R 
Eq. (E2D, with 7 2 = 4ft A . 



and Qr = 
From (J23D 

(25) 



6(2H 2 + iJ) from the above expressions yields 



In Appendix A we provide a detailed analysis of the vacuum MC equations with 
Qm = 0. Although the presence of matter plays an important role in our minimal 
curvature construction, an analysis of the vacuum theory provides valuable insight into 



the solutions we are interested in studying. In |Appendix B[ we transform the vacuum 
MC theory into an Einstein frame and relate quantities of physical interest in both 
frames. 

For general solutions to the equations of motion ( |20|) and ( 1211 with functions ( |T7I) 
we solve for the Hubble parameter H(t) and scalar field (p(t). We plot the relevant 
portion of the (p—H phase space in Fig. [TJ 



2 . 5 



2.25- 



1 .75 



1.5- 



1.25 - 




0.5 



1.5 



2.5 



Figure 1. Phase space of ip vs. H . All solutions asymptotically approach a late-time 
accelerating phase with constant H at the minimal curvature 1Z denoted by the red 
line. 

To solve the equation we integrate from the past to today (H = 1) and then 
from today into the future and patch the solutions together. In the plots we take the 
conditions H Q = 1.0, 7 = 3.15, Qm = 0.25 and let the values of ipo vary. At late times, 
the solutions approach the constant H attractor when the Universe is well into the 
accelerating epoch. 
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3. Comparison with ACDM 

We now compare our model with ACDM. To do so, we must enter reasonable initial 
conditions into our numerical study and solve the equations of motion fl20l) . fl2~Tl) together 
with ( fl7l) and the equation a = Ha. Let us begin by considering the value of the minimal 
curvature. Physically the minimal curvature corresponds to an emptying of the matter 
in the Universe due to cosmological expansion. In our model the value of the minimal 
curvature is approached asymptotically. Today, the value of the curvature is given by 

i?o = 6(2# 2 + if )~12# 2 . (26) 

Because we observe a significant amount a matter in the Universe today, we know that 
we have not yet reached the minimal curvature. However, because we are accelerating, 
we know that we are near the minimal curvature (i.e. the first and second terms on the 
right hand side of equation ( 1221) must be comparable). Hence, from equation ( 1261) . we 
expect the value of the minimal curvature 1Z to be close to but less than 1Z < 12Hq. 
Therefore, in terms of our dimensionless quantities ( fl8l) . 1Z < 12 and the free parameter 
in our theory 7 < \[V1 fa 3.4641. For the solutions considered below, we choose a 
value of 7 = 3.15 meeting the above requirements and that follows the ACDM model 
to our satisfaction for our toy construction (i.e. a value leading to a matter dominated 
cosmology followed by a "jerk" near a redshift of z ~ .5 into an accelerating phase). 
In the action (jBJ), the effective Newton's constant is GN eff is given by 

{lftoG NtJJ )-i = -±- rp . (27) 

To ensure that the effective Newton's constant remains positive definite over the history 
of the Universe ( GN eff > 0, Vt) we must have 7^ < .5. We are almost in a position to 
compare our model both with ACDM, which has the Friedmann equation 



H 



\ 



(o) 

M , m n (°) 



a 3 



+ (i-nff), (2* 



and with the observational data provided by the SNela gold sample. To make a 
comparison with the observational data we require an understanding of the luminosity 
distance formula in the context of modified gravity models. 

An important consideration arises when using the formula for the luminosity 
distance in theories of the form: 



S 



J d A x^g 



(29) 



where the £j represent the different types of possible matter Lagrangians present. Such 
a theory arises as the low-energy effective action for the massless modes of dilaton 
gravity in string theory, and our model is an action of just this sort [H1H2]; albeit, with 
an unusual choice for the functions A, B, V, C . As we have already discussed, these 
theories typically lead to time variation in Newton's gravitational constant. The time 
variation can affect the way one should compare the theory to observations [50J. In 
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particular, the time-evolution can alter the basic physics of supernovae. For example, 
the time variation in Gn leads to different values of the Chandrasekhar mass at 
different epochs, and hence, a supernova's peak luminosity will vary depending on when 
the supernova exploded. This makes treating the supernovae la as standard candles 
difficult [5TJ E2, E31 EU EE, EE] . Specifically, the peak luminosity of SNela is proportional 
to the mass of nickel synthesized which is a fixed fraction of the Chandrasekhar mass 
Mch ~ G~ 3 / 2 . Hence, the luminosity peak of SNela varies as L ~ G~ 3 ^ 2 and the 
corresponding absolute magnitude evolves as 

M = M + ^log^, (30) 

4 Cr 

where the subscript zero indicates the local values of the quantities. Therefore, the 
magnitude-redshift relation of SNela in modified gravity theories of the type given by 
( |29l) is related to the luminosity distance via [53, 55J: 

m(z) = M + 5 log d L {z) + H log . (31) 

4 Cr 

Even if gravitational physics is described by some theory other than General Relativity 
the standard formula for the luminosity distance applies as long as one is considering a 
metric theory of gravity |57j : 



^ (32) 

For ACDM, the Luminosity distance (1321) can be written [40J: 
dh(z\ ^a, Hq) = 

' r [l + z') 2 {l + Q M z')-z'{2 + z')Q A \^) (33) 



c(l + z) 

S \ \k\ 



Hn\ \k\ V ^ 



where, S(x) = sin(x) and k — 1 — Vt tot for Q tot > 1 while S(x) = sinh(x) with 
k = 1 — Q t ot for fltot < 1 while S(x) = x and k = 1, for Q tot = 1. Here and throughout, 

^tot = + ^A- 



3.1. Numerical Analysis 

Given the above considerations we take the following values of parameters today: 

a = 1.0, H = 1.0, cp = 0.1, 7 = 3.15, = 0.25 . (34) 

We then integrate our equations from the past to today and then from today into the 
future and patch the solutions together. While we do not provide an exhaustive study of 
the parameter space for the MC model, the parameters given above provide a successful 
example of the construction which fulfills our rather modest goals 0. It is quite possible 
that the parameters (I34"]) can be tuned to achieve an even better agreement with ACDM. 

* An exhaustive study of the phase space of the vacuum theory is supplied in |Appcndix A| 
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The history and future of the curvature R together with the co-moving Hubble radius 
are plotted in Fig. El 




Figure 2. The Ricci scalar (left) and co-moving Hubble radius (right) as functions of 
cosmic time t. The curvature decreases from the matter-dominate past (red curve) to 
the constant minimal value at R = j 2 indicated by the green line. Today, to = 1 and 
the blue curve indicates the future evolution of the curvature. While the Universe is 
accelerating the co-moving Hubble radius decreases. 



In the Figure the past history of the curvature is plotted in red while the future is 
plotted in blue. Today we are at the value t — 1. The minimal curvature is given by 
the green line at R = TZ = 7 2 . As expected, the curvature is large in the past when 
the Universe is matter dominated and then decreases, approaching the accelerating 
late-time de-Sitter attractor with constant minimal curvature R = 1Z. The second 
plot in the Figure shows the evolution of the co-moving Hubble radius H~ 1 /a. The 
phenomenologically desired cosmological transition from matter domination to late- 
time acceleration of the Universe is clearly indicated by the decreasing of the co-moving 
Hubble radius, when 

~ <0. (35) 
at a 

In some modified gravity models it can be difficult to achieve this transition (see, 
e.g. [5EJE9]). 

The majority of our results are presented in Figures [3] and HI In Figure [3] we plot 
the scale factor a as a function of cosmic time t, the Hubble parameter H as function 
of redshift z, where 

z + l=^, (36) 
a 

the deceleration parameter 

«M = -Ji. (37) 
and the luminosity distance, d^. 
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a.) Scale Factor _ b.) Hubble Parameter 



q( 




. 2 50 . 50 . 75 1 1.251.5 



z 

Figure 3. A comparison of the cosmological evolution of the MC model with ACDM. 
Quantities plotted in red are the MC model while those in blue are for ACDM. Today 
we are at the values t = 1, z = 0. Plot a.) shows the time evolution of the scale factor. 
Plot b.) shows the Hubble parameter as a function of rcdshift. Plot c.) compares 
the deceleration parameters as a function of redshift. Plot d.) shows the luminosity 
distance, d^, as a function of redshift. The MC model is clearly a good fit with ACDM. 

For the particular set of parameters considered, the scale factor of our model 
differs from ACDM most strongly in the far future, although the difference in the 
expansion rates of the MC construction with ACDM is apparent in the plot of the Hubble 
parameters at high redshifts. The transition from matter domination to acceleration (the 
jerk) occurs at a slightly lower redshift than ACDM. There is only a slight difference 
in di that occurs at high redshifts (z > 1), although the difference is not significant 
enough to distinguish our model from ACDM using only the current supernova data. 
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Figure 4. Comparison of the MC and various ACDM models with observational data. 
The supernovae data points are plotted with error bars and the data is taken from [27j . 
The luminosity distance d^ for the MC model is plotted by the dashed (blue) curve. 
The various theoretical predictions for ACDM are represented by the solid curves and 
were examined in 281. 



In Figure HI we compare the luminosity distance predicted by our model with several 
versions of ACDM and with the observational data. The luminosity distances are plotted 
out to a redshift of z = 2 (the highest redshift supernova data is from z = 1.76). The 
theoretical predictions of the minimal curvature construction and ACDM are compared 
with the "gold" supernovae sample of [27]. The particular choice of ACDM models 
shown are from [28]. The luminosity distance for the minimal curvature construction 
is denoted by the blue dashed line and fits the supernova data extremely well. The 
luminosity distance cIl is virtually indistinguishable from ACDM with = 0.25 and 
O4 0) = 0.75 (there is a small discrepancy at high redshift as can be seen from Figure [3]). 



4. Conclusions 



We have shown that a period of late-time cosmic acceleration can follow directly 
from a simple minimal curvature conjecture (MCC). The model fits the SNela data 
exceptionally well. While the specific formulation considered here is only a toy 
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construction, unlikely to be compatible with constraints from solar system and table 
top test of the equivalence principle, it may capture phenomenologically interesting 
features of a more fundamental theory that admits a limiting minimal curvature. 
Furthermore, the construction successfully demonstrates the possibility that a new- 
fundamental physical principle may ultimately be responsible for the recent period of 
cosmological acceleration. It is possible that experimentally viable models based on the 
minimal curvature conjecture exist. The search for such models within the context of 
scalar-Gauss-Bonnet gravity is currently underway. Despite the tight theoretical and 
experimental constraints on scalar- Gauss- Bonnet cosmologies [62J E3 Ell ESI ES], we 
remain optimistic that an experimentally and theoretically viable model based on the 
minimal curvature construction can be discovered. 



Acknowledgments 

It is a pleasure to thank R. Brandenberger, R. Gregory, V. Jejjala, I. Moss, R. Myers 
and T. Underwood for helpful discussions. I am especially grateful to M. Trodden for 
numerous useful discussions over the course of this work. This work is supported in part 
by PPARC and by the EU 6th Framework Marie Curie Research and Training network 
"UniverseNet" (MRTN-CT-2006-035863) . 



References 



A. G. Ricss et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998); 

[arXiv:astro-ph/9805201| ; 
S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999); 

[ar Xiv : astr o-ph/ 9 8 1 2 1 3 3| ; 
C. B. Netterfield et al. [Boomerang Collaboration], Astrophys. J. 571, 604 (2002); 



arXiv: astro-ph/0 1 04460 : 
N. W. Halverson et al., Astrophys. J. 568, 38 (2002) . |arXiv:astro-ph/0104489l . 
J. L. Tonry et al., Astrophys. J. 594, 1 (2003); |arXiv:ast ro-ph/03050 08[r~ 
C. L. Bennett et al., Astroph ys. J. Su ppl. 148, 1 (2003) ; |arXiv:astro-ph/0302207] ; 
S. Nojiri and S. D. Odintsov, |arXiv:hep-th/0601213| 

E. J . Copeland, M. S ami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) 
[arXiv:hep-th/0603057| . 

S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, Phys. Rev. D 70, 043528 (2004) 



arXiv:astro-ph/0306438|. 



S. M. Carroll, A. De Felice, V. Duv vuri, D. A. E asson, M. Trodden and M. S. Turner, Phys. Rev. 
D 71, 063513 (2005) [arXiv:astro-ph/041003l| . 

D. A. Easson, Int. J. Mod. Phys. A 19, 5343 (2004) [arXiv:astro-ph/0411209| . 

S. Nojiri and S. D. Odintsov, Mod. Phys. Lett. A 19, 627 (2004) |arXiv:hep- th/03100 45). 
F. P. Schuller and M. N. R. Wohlfarth, Phys. Lett. B 612, 93 (2005) arXiv:gr- qc/0411076|. 
F. P. Schuller and M. N. R. Wohlfarth, Nucl. Phys. B 698, 319 (2004) ar Xiv:hep-t h/0403056]. 
M. A. Markov, Pisma Zh. Eksp. Teor. Fiz. 36, 214 (1982); Pisma Zh. Eksp. Teor. Fiz. 46, 342 

(1987); V. P. Frolov, M. A. Markov and V. F. Mukhanov, Phys. Lett. B 216, 272 (1989). 
V. P. Frolov, M. A. Markov and V. F. Mukhanov, Phys. Rev. D 41, 383 (1990). 
D. Morgan, Phys. Rev. D 43, 3144 (1991). 

V. Mukhanov and R. H. Brandenberger, Phys. Rev. Lett. 68, 1969 (1992); 

M. Trodden, V. F. Mu khanov and R. H. Brandenberger, Phys. Lett. B 316, 483 (1993); 
|arXiv:hep-th/9305 111] ; 

[20] R H. Brandenberger V. Mukhanov and A. Sornborger, Phys. Rev. D 48, 1629 (1993); 
|arXiv:gr-qc/9303001] . 



The Accelerating Universe and a Limiting Curvature Proposal 15 



R. Moessner and M. Trodden, Phys. Rev. D 51, 2801 (1995); [arXiv:gr-qc/9405004]. 

R. H. Brandenberger, R. Easther and J. Maia, JHEP 9808, 007 (1998); arXiv:gr-qc/9806111 . 
D. A. Easson and R. H. Brandenberger, JHEP 9909, 003 (1999); |arXiv:hep-th/9905175|. 
D. A. Easson and R. H. Brandenberger, JHEP 0106, 024 (2001); |arXiv:hep-th/0103019|. 
D. A. Easson, JHEP 0302, 037 (2003); |arXiv:hep-th/0210016l; 
D. A. Easson, Phys. Rev. D 68, 043514 (2003). |arXiv:hep-th/0304168| . 

A. G. Riess et nl. [Supernova Search Team Collaboration], Astrophys. J. 607, GG5 (2004) 

[arXiv:astro-ph/0402512| . 
T R. Choudhury and T. Padmanabhan, Astron. Astrophys. 429, 807 (2005) 

[arXiv:astro-ph/0311622| . 
G. Magnano and L. M. Sokolowski, Phys. Rev. D 50, 5039 (1994) [arXiv:gr~qc/9312 008j. 
T. Damour, G. W. Gibbons and C. Gundlach, Phys. Rev. Lett . 64, 123 (1990). 
R. Bean and J. Magueijo, Phys. Lett. B 5 17, 177 (2001) |arXiv:astro-ph/0007199| . 
R. Bean, Phys. Rev. D 64, 12 3516 (2001) |arXiv:astro-ph/0104464| . 
T. Biswas and A. Mazumdar, | ar Xiv:hep- ty 408026. 

T. Bi swas, R. Brandenbe rger, A. Mazumdar and T. Multamaki, Phys. Rev. D 74, 063501 (2006) 



arXiv:hep-th/ 0507199|. 



R. Bean, S. Carroll and M. Trodden, arXiv:astro-ph/0510059. 
S. Das, P. S. Corasaniti and J. Khoury, Phys Rev. D 73, 083509 (2006) [arXiv:astro-ph/0510628] . 
L. Amendola, D. Polarski and S. Tsujikawa, |arXiv:astro-ph/0603703. 



S. M. Carroll, I. Sawicki, A. Silvestri and M. Trodden, arXiv:astro-ph/0607458 
J. P. Uzan, Rev. Mod. Phys. 75, 403 (2003) [arXiv:hep^ph /020534b"[7^ 

S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 483, 565 (1997) 

[arXiv:astro-ph/9608192] . 
T. Damour and A. M. Polyakov, Gen. Rel. Grav. 26, 1171 (1994 1 |arXiv:gr-qc/ 9411069j. 
T. Damour and A. M. Polyakov, Nucl. Phys. B 423, 532 (1994) |arXiv:hep-th/9401069|. 
D. A. Easson, Int. J. Mod. Phys. A 16, 4803 (2001) |arXiv:hep -th/0003086 . 

D. Langlois, Prog. Theor. Phys. Suppl. 148, 181 (2003) |arXiv:hep-th/0209261| . 

P. Binetruy, C. Deffaye t, U. Ellwanger and D. Langlois, Phys. Lett. B 477, 285 (2000) 

[arXiv:hep-th/9910219| . 
Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (200 6) |arXiv:astro-ph/06040511 . 
M. Doran and G. Robbers, JCAP 0606, 026 (2006) |arXiv:astr o-ph/0601544p" 

E. V. Linder, Astropart. Phys. 26, 16 (2006) arXiv:astro-ph/0603584 . 

D. Langlois, Prog. Theor. Phys. Suppl. 163, 258 (2006) |arXiv:hep-th/0509231| . 

V. Acqua yiva, C. B accigal upi, S. M . Leach, A. R. Liddle and F. Perrotta, Phys. Rev. D 71, 104025 

(2005) [arXiv:astro-ph/0412052] . 
L. Amendola, P. S. Corasaniti and F. Occhionero, arXiv:astro-ph/9907222. 



E. Garcia-Berro, E. Gaztanaga, J. Isern, O. Benvenuto and L. Althaus, arXiv:astro-ph/9907440 
A. Riazuelo and J. P. Uzan, Phys. Rev. D 66, 023525 (2002) |arXiv:astro-ph/0107386p 
E. Gaztanaga, E. Garcia-Berro, J. Isern, E. Bravo and I. Dominguez, Phys. Rev. D 65, 023506 

(2002) [arXiv:astro-ph/0109299| . 
S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 73, 103511 (2006) [arXiv:astro-ph/0602053| . 
R. Gannouji, D. Po larski, A. Ranquet and A. A. Starobinsky, JCAP 0609, 016 (2006) 
arXiv:astro-ph/0606287| . 



C. Shapiro and M. S. Turner, |arXiv:astro-ph/0512 586. 

D. A. Easson, F. P. Schuller , M. Trodden and M. N. R. Wohlfarth, Phys. Rev. D 72, 043504 (2005) 
[arXiv:astro-ph/0506392| . 

R. Punzi, F. P. Schuller and M. N. R. Wohlfarth, arXiv:gr-qc/0605017[ 

H. K. Jassal, J. S. Ba gla and T. Padmanabhan, Phys. Rev. D 72, 103503 (2005) 

[arXiv:astro-ph/0506748| . 

H. K. Jassal, J. S. Bagla and T. Padmanabhan, [arXiv:astro^ ph/0601389[ 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 631, 1 (2005) ~ [aildv:hep-th/0508049] . 

G. Cognola, E. Elizalde, S . Nojiri, S. D. Odintsov and S. Zerbini, Phys. Rev. D 73, 084007 (2006) 

[arXiv:hep-th/0601008| . 

A. De Felice, M. Hindmarsh and M. Trodden, JCAP 0608, 005 (2006) |arXiv:astro-ph/0604154] 
G. Calcagni, B. de Carlos and A. De Felice, Nucl. Phys. B 752, 404 (2006) |arXiv:hep-th/0604201|. 
T. Koivisto and D. F. Mota, Phys. Lett. B 644, 104 (2007) |arXiv:astro -ph/0606078 . 



The Accelerating Universe and a Limiting Curvature Proposal 



16 



Appendix A. Probing the Vacuum Theory 

The basic mechanism used to implement the limiting curvature construction requires a 
non-minimal coupling of (p to the matter Lagrangian. This coupling plays an important 
role in the overall cosmological dynamics discussed above. Hence, an understanding of 
the vacuum theory is perhaps not as fruitful as it would be in case of a theory with 
minimal coupling of ip to matter; however, once the Universe becomes sufficiently large 
and dilute the dynamics will resemble those of the vacuum theory As we shall see the 
vacuum MC theory is significantly richer than that of ACDM. The equations describing 
the vacuum are given by the equations ( fl6|) with Qm = 0. The dimensionless forms are 

H 2 - 2-ipH 2 - 2-jHp + = o , 

3 

7 (2tf 2 + if)-^M = 0, (A.l) 
o 

where we have re-introduced the parameter \i from equation (|T7|) (previously set equal 
to 7) for completeness and the potential is V(<p) = fi 3 <p. In the case of the vacuum we 
may solve exactly for H(t). We focus on two types of solutions of particular interest. 
In the first, the Hubble parameter is given by 



H(t) = tanh - 6a 7 ) } , (A.2) 

where72 = /i 3 /7 is the minimal curvature and a is a constant. In this case 

?(t-6a 7 )>, (A.3) 




and the scale factor evolves as 
a(t) = a 




soeh { \/^-(6n - — t) 



(A.4) 



The quantities mentioned above, along with the field <p are plotted in Fig. (1A1I) . In this 
case the Universe undergoes a bounce as indicated in plots a.) and b.) in Fig. ( 1A1I) . 
After the bounce the Hubble parameter is pulled up to the value at the minimum 



curvature H — > JlZ/12. This is a de Sitter solution of the vacuum theory with 



H dS = JTZ/12, (A.5) 



while cp continues to evolve exponentially <£>(i) + (27) 1 oc eaq>(y 72-/12). Interestingly, H 
and H evolve in such a way that the curvature R = 6(2H + H) is constant throughout 
the evolution of the solution, fixed at the minimal value R = 1Z. This behavior is not 
surprising as the constraint on R comes from the second equation in the EOM (IA.1I) . 

The second family of solutions are of greater relevance. They are the vacuum 
analogs of the solutions plotted in Fig. ^ and discussed at the end of Section [2J These 
solution are pulled down to the value of H at the minimal curvature. The entire phase 
space of solutions to the vacuum theory is shown in Fig. flA20 . Representative solutions 
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Figure Al. Various quantities of interested plotted as functions of cosmic time t. In 
the plots we take lZ = j = (i = a = l. 



of the two families of solutions discussed above are plotted by the solid curves. The 
red curve is the solution given by equation ( 1A.2I) . The second set of solutions (relevant 
to Section [2]) are plotted in green in the upper right quadrant of the phase space. 
The system defined by the vacuum equations (equation ( lA.lj) ) has two unstable saddle 
equilibrium points at the values 




(A.6) 



The equilibrium points are marked by the purple dots in Fig. (1A2I) . 
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Figure A2. The phase space of the vacuum theory. H is plotted on the vertical axis 
and if is plotted on the horizontal axis. The vector field for solutions is drawn in black. 
Solutions of particular interests are plotted by the solid curves. The constant de Sitter 
value of H at the minimal curvature is indicated by the dashed (blue) line. The two 
unstable saddle equilibrium points are designated by the purple dots. In the plot we 
take 7 = /i = 3.15. 



Appendix B. An Einstein Frame Description 

It is not our intent to provide a complete Einstein frame analysis of the full MC system 
with matter; however, it is instructive to consider the vacuum theory in an Einstein 
frame. To move to the Einstein frame we begin with the Brans-Dicke frame defined by 
the action ([8]) . Passage to the Einstein frame is achieved via a conformal transformation 
of the form 

g^ = VL 2 g^ V} (B.l) 

where Q 2 is the conformal factor which must be positive to leave the signature of the 
metric unaltered. From this point forward a tilde shall denote a quantity built out of 
the Einstein- frame metric tensor g^ v . Under this transformation the infinitesimal line 
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element and the determinant of the metric transform as 

ds 2 = tt 2 ds 2 , J~-g = tt A y/=g~ , (B.2) 
respectively. The Ricci scalar transforms as 

R = Q 2 (r+ 6D(lnn) - 6^ V ^J^ . (B.3) 

Omitting the ordinary divergence y/— pQ(ln Q), the action in the Brans-Dicke frame 
transforms to 

S = \d\f^g (JLr - e^V^fi + , (B.4) 

where the potential in terms of <3> is given by ([7j) together with f fT7|) : 

V($) = ^^-$j . (B.5) 

Here we have re-introduced the parameter \i from equation ( TT7T) (previously set equal 
to 7) for completeness. 

By choosing our conformal factor to be Q 2 = 2$/M 2 1 , and performing a field 
redefinition to define the Einstein frame field d>, 



*=^ex P « tJL. , (B.6) 




the action (1B.4h becomes the Einstein frame action 



Sef = Jd A x^g \ ~~2~R ~ 2^V M 0V,0 - V{<j>) J , (B.7) 
where we have defined the potential 



$ 

The Einstein frame potential is given by 

M 2 

V (0) = ^Bift (e w ^' - l) e" 2 ^^ , (B.9) 

where 7?. = /i 3 /7 is the minimal curvature and (3 = v/2/3. 

In this frame the field (/> has a canonically normalized kinetic term making the 
interpretation of solutions more simple due to our familiarity with minimally coupled 
scalar field to ordinary Einstein gravity. It is important to note, however, that so far 
we have only considered the vacuum theory. Even in the Einstein frame the Einstein 
field will couple non-minimally to matter and therefore, when matter is present in 
significant amounts, the simple Einstein frame vacuum solutions will not be an accurate 
description of the theory. 

Under the conformal transformation ( IB.lj) . the cosmic time coordinate transforms 

as 

dP = e^ M ^dt 2 . (B.10) 
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Taking the Einstein-frame Friedmann-Robertson- Walker flat metric 

ds 2 = -dt 2 + d 2 (t)dy? , 
Leads to the familiar equations of motion 
H 2 



and 



fr" + 3H<p' + 



dV 
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(B.H) 
(B.12) 

(B.13) 



where prime denotes differentiation with respect to the Einstein frame cosmic time 
coordinate t and H = a! fa. Using the conformal transformation fIB.lj) . we find the 
Hubble parameters in the Einstein and Minimal Curvature frames are related by 



H 



P<t> 



'pi H 







2M, 



pi 



(B.14) 



We plot the Einstein frame potential flB.9j) in Fig. fIBll) . From the plot it is clear 
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Figure Bl. The Einstein frame potential V as a function of the Einstein frame field 
4>. In the plot we have set WttGn = 7 = /i = 1. 

that there are three regions of interest: 

i) There is an unstable de Sitter solution corresponding to the field sitting on the 
top of the potential. 

ii) If the field starts initially to the left of the unstable de Sitter solution (and 
without significant positive velocity) the system will quickly run off to large negative <fi. 

iii) If the field falls to the right of the de Sitter solution it will asymptotically roll 
to large positive values. 
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Let us examine the first possibility in greater detail. The unstable de Sitter solution 
with constant H = H dS sits at the value at the maximum of the potential (IB. 91) . The 
maximum is at constant value 

</> = <f> dS = p- 1 M&]n2, (B.15) 

From equation (IB. 121) we find the corresponding H 

H*s=lM- (B.16) 



From equation (1B.14|) we see that the corresponding value in the minimal curvature 
frame is 

H(H dS ,cf> dS ) = ^, (B.17) 

which is the exact de Sitter solution we found in our analysis of the vacuum of the 
minimal curvature theory (|A.5jl . 

The Minimal Curvature frame field is related to the Einstein frame field via 

M 2 

^)3(l- e ^), (B.18) 
and consequently, 

* w = X b, ( 1 -^)- (R19) 

Using the relation (1B.18I) . we find the value of (p at the de Sitter point (IB. 151) : 

Vttds) = ■ (B.20) 

Hence we conclude that the unstable de Sitter solution in the Einstein frame is mapped 
to one of the unstable saddle critical points in the vacuum of the Minimal Curvature 
frame (see equation flA.6j) ) . 

We now examine case iii): when rolls to large positive values the Einstein frame 
potential may be approximated by 

M 2 7? 

V(d>) ~ iVV/W , (B.21) 
leading to the exact solutions 

m -VQ\n{\l^—\ (B.22) 




M pl 

H(i) = I , (B.23) 

where Vq = M^{R/2, corresponding to power- law acceleration in the Einstein frame with 

a(t) = a P . (B.24) 

To find the corresponding behavior in the Minimimal Curvature frame we use the 
conformal transformation (IB. ID along with the transformation for the cosmic time 
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coordinate (IB.lOp . We find the late-time power law attractor in the Einstein frame 
( IB. 241) is mapped to the asymptotic de Sitter attractor in the MC frame at the minimal 
curvature 71 (depicted by the dashed blue line in Fig. ( 1A2I) ). 

a(t)=a e^ t . (B.25) 



